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Abstract:  Precise  calculations  of  absorbed  dose  (AD)  are  a difficult  task  involving  physical  phenomena  such  as  emission, 
transport,  and  absorption  of  radiation.  We  used  the  Monte  Carlo  method  to  calculate  Kerma  as  well  as  AD  in  water  imparted  by 
two  different  lr-192  HDR  brachytherapy  seeds  ( Flexisource  and  microSelectron-v2)  taking  into  account  the  AAPM  TG-43 
formalism  and  AAPM  & ESTRO  most  recent  reports  recommendations.  The  aim  of  this  work  is  to  evaluate  when  can  Kerma 
be  used  as  a measurement  of  AD  for  this  type  of  seeds.  Thus,  we  analyse  the  behaviour  of  both  quantities  in  whole  space, 
putting  special  emphasis  near  the  source  surface.  We  carried  out  calculations  using  microvoxels  to  obtain  high  spatial 
resolution  of  data  close  to  the  source.  We  observed  differences  of  up  to  6 % between  AD  and  Kerma  within  1 mm  around  the 
seed,  and  less  than  1 %>  in  any  other  region  of  the  phantom.  This  allows  us  to  analyse  build-up  region  for  lr-192  HDR 
brachytherapy  seeds.  As  it  will  be  further  discussed  in  this  paper,  our  results  can  be  explained  in  terms  of  partial  electronic 
equilibrium  reached  on  different  regions  of  the  phantom.  Both  seeds  showed  common  overall  behaviour,  providing  generality 
to  the  conclusions  drawn.  The  complete  bearing  of  the  radial  dose  function  (defined  in  the  TG43  formalism)  as  it  traverses  the 
surface  of  the  seed  is  reported.  Whenever  comparisons  are  possible,  our  results  are  in  agreement  with  those  reported  by  other 
authors.  Tables  of  radial  dose  function,  including  new  data  computed  from  AD  rate  (instead  of  Kerma  rate),  are  presented. 
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1.  Introduction 

In  HDR-brachytherapy  treatment,  the  radioisotope  lr-192 
is  used  due  to  its  specific-activity  and  the  characteristics  of 
its  photon-spectrum  in  a sealed  source  form.  For  use  in 
treatment,  these  sources  are  manufactured  in  miniaturized 
size  of  approximately  cylindrical  shape  (about  1 mm  diameter 
by  several  millimetres  long,  see  Fig.  1)  and  they  are  known 
as  seeds.  The  seeds  may  vary  slightly  in  design  depending  on 
the  manufacturer  [1].  All  HDR-brachytherapy  seeds  have  a 
steel-shielding  that  covers  their  iridium  core.  This  shielding 
provides  rigidity  to  seeds,  being  designed  to  stop  beta- 
particles  from  radioactive  decay  while  allowing  the  outward 
of  photons  flux. 

During  treatment,  the  seed  is  inserted  into  the  patient  using 
an  applicator.  Radiation  emitted  from  the  source  (i.e.  primary 
photons)  interacts  with  the  surrounding  environment  (human 
tissue)  setting  in  motion  secondary  electrons.  Along  their 
path,  the  kinetic  energy  of  the  secondary  electrons  is 
absorbed  by  the  medium.  By  means  of  this  process, 
radiation-dose  is  absorbed  by  the  tissues. 


In  regions  of  space  where  electronic  equilibrium  is  attained, 
AD  (absorbed  dose)  from  primary-photons  can  be  inferred 
from  knowledge  of  Kerma  (kinetic  energy  released  per  unit 
mass).  In  contrast,  near  seed  surface  there  is  a build-up 
region  where  electronic  equilibrium  is  not  achieved  and 
marked  differences  between  Kerma  and  AD  are  expected. 
Other  authors  [4-6]  have  attempted  to  characterize  this  build- 
up region  in  the  past.  Although,  as  it  will  be  further  discussed 
in  this  paper,  calculations  achieved  by  them  are  based  in 
models  of  simplified  situations.  We  have  not  found  precise 
reports  of  direct-computation  of  AD  over  the  entire  region 
around  the  seeds.  The  use  of  Kerma  as  a measure  of  AD  from 
photons  remains  a current  topic  of  discussion. 

In  addition,  precise  knowledge  of  the  AD  distribution  in 
water  is  required  in  order  to  plan  the  treatment.  The  Monte 
Carlo  (MC)  simulation  method  showed  being  an  appropriate 
tool  for  calculating  AD  distribution  in  a medium  around  the 
seed  [7].  The  MC  implementation  to  this  problem  requires 
following  a set  of  well-established  protocols  [1,  8-9]. 
Calculating  AD  is  more  complex  than  Kerma  computation 
and  consequently,  it  requires  increased  computational  time. 
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Thus,  it  is  important  to  determine  precisely  in  which 
situations  AD  can  be  inferred  from  Kerma  computation,  this 
being  one  of  the  goals  of  this  work. 

It  is  important  to  state  that  in  the  past,  many  authors  have 
applied  the  MC  method  to  model  HDR-brachytherapy  seeds 
with  several  degrees  of  approximation  and  to  calculate 
Kerma  in  a water-phantom:  microSelectron-v2  [3,  10-11], 
microSelectron-v2r  [12]  (slightly  different  from  v2  seed)  and 
Flexisource  [2,  11].  According  to  the  literature  on  this  subject, 
to  calculate  AD  even  closer  to  the  seeds,  authors  frequently 
have  resorted  to  data  extrapolation  [1]. 

2.  Materials  and  Methods 

In  this  work  we  used  the  MC  PENELOPE  v2008  code  [13- 
14]  together  with  the  clonEasy  simulation  package  available 
at  http://inte.upc.edu/downloads  [15].  Two  available 
commercial  lr-192  HDR-brachytherapy  seeds  have  been 
modelled:  Flexisource  and  microSelectron-v2.  Radiation 
transport  in  water  has  been  simulated  following 

recommendations  given  by  AAPM  in  TG43  formalism  [8] 
and  its  update:  TG43U1  [9]  and  new  AAPM  and  ESTRO 
report  [1].  The  photons  interactions  considered  in  the 


PENELOPE  MC  code  are:  Rayleigh  dispersion,  photoelectric 
absorption,  Compton  dispersion  and  pair  production.  The 
photoelectric  absorption  cross-sections  used  by  PENELOPE 
are  interpolated  from  the  LLNL  evaluated  photon  data  library 

[16]  (EPDL).  Differential  cross-sections  (DCS)  per  unit  solid 
angle  for  Rayleigh  dispersion  are  calculated  using  non- 
relativistic  perturbation  theory  with  atomic  form  factors 
obtained  from  EPDL  library  [16].  Compton  dispersion  is 
described  according  to  the  relativistic  impulse  approximation 

[17]  where  the  Doppler  broadening  effect  is  taken  into 
account.  For  pair-production  a semi-empirical  model  based  in 
Bethe-Heitler  DCS  equation  is  used  [18].  As  regards 
electrons,  PENELOPE  implements  a combined  detailed- 
simulation  of  hard  events  and  condensed-simulation  of  soft 
events  in  a mixed  scheme  [19-21].  Concerning  electrons  and 
positrons,  elastic  scattering  cross-sections  are  computed 
using  interpolation  of  numerical  DCS  from  the  ELSEPA 
database.  The  simulation  of  inelastic  collision  of  electrons 
and  positrons  is  performed  according  to  Bethe  and  Fano  DCS 
equation,  and  on  the  basis  of  an  atomic  generalized  oscillator 
model.  The  DCS  for  bremsstrahlung  emission  is  based  on  the 
Seltzer  and  Berger  original  database  [19,  20]. 


Figure  1.  Modelled  HDR  brachytherapy  Ir-192  seeds:  (a)  Flexisource  and  (b)  microSelectron-v2.  Geometry  dimensions  (in  mm)  and  composition  of 
commercial  seeds  were  taken  from  literature  (Flexisource  [2],  microSelectron-v2  [3]).  Simulated  homogeneous  nuclear  distribution  for  each  of  the  seeds 
modelled  in  this  work  is  also  included  (in  low  statistic,  only  for  visualizing  purposes). 


Both,  seeds  composition  and  geometry  have  been  taken 
from  the  literature:  Daskalov  et  al  [3]  for  microSelectron-v2 
seed  and  Granero  et  al  [2]  for  Flexisource  seed.  Active  core 
was  modelled  by  a homogeneous  iridium  atoms  distribution 
taking  into  account  precise  cylindrical  geometry  of  each  one 
of  the  seeds  (see  Fig.  1).  The  seeds  were  modelled  including 
the  stainless  steel  capsule  and  flexible  cable  for  remote 
positioning.  The  3-D  space  in  which  all  physical  processes 
were  simulated  consisted  of  a spherical  water  phantom  (40 


cm  radius)  with  the  seeds  positioned  at  its  centre. 
Discretization  of  the  region  where  Kerma  and  AD  are 
calculated  is  achieved  into  cubic  voxels  with  varying  edge- 
lengths:  0.01  mm,  0.1  mm  or  1 mm  (increasing  with  seeds 
surface  distance).  This  high  spatial  resolution  compromises 
the  minimum  number  of  simulated  histories  necessary  to 
achieve  statistically  representative  data.  In  addition,  studying 
dose  high  spatial-gradient  without  introducing  additional 
error  requires  considering  voxel  dimensions  appropriated  to 
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this  end.  Otherwise,  calculation  of  the  near-seed  dose  could 
be  over-  or  under-estimated  in  calculations.  Table  1 displays 
data  from  the  literature  [2-3]  that  has  been  used  in  this  work 
to  model  seeds  and  phantoms. 


Table  I.  Main  geometry  dimensions  of  phantoms  and  seeds,  and  densities 
and  composition  by  weight  of  each  material  simulated  in  this  work:  (a) 
Flexisource  [2]  and  (b)  microSelectron-v2  [3]. 


Figure  2.  Energy  spectrum  corresponding  to  photons  emitted  by  Ir-192  in 
100  decayment  events  to  Os -192  (EC,  5 %)  and  Pt*-192  (fd  ~,  95  %)  isobars 
in  nuclear  and  atomic  excited  states.  In  embedded  window,  photon  groups  of 
greater  intensity  are  appreciated. 


In  order  to  use  the  TG-43  formalism  and  compare  Kerma 
with  photons-AD,  we  have  considered  only  X-  and  y-rays  as 
primary  radiation  [24],  without  including  (3  particles,  Auger 
or  CE  (conversion  electrons,  see  Fig.  2).  Such  particles  are 
largely  blocked  by  the  steel-shielding  that  constitutes  the 
capsule  of  the  seeds  (see  Fig.  1).  Besides  this,  the  main  goal 
of  this  work  is  to  evaluate  validity  of  Kerma  as  a measure  of 
water  AD  from  photons.  In  AD  computation,  simulation 
process  of  each  history  ends  when  average  kinetic  energy  of 
all  set-in-motion  electrons  ensure  this  energy  is  deposited  in 
a voxel,  and  represents  cut-off  energy.  We  set  cut-off  energy 
of  103  eV  in  this  work  simulation.  This  cut-off  energy  value 
ensures  that  the  linear  dimensions  of  the  voxel  are  greater  by 
an  order  of  magnitude  than  the  mean-range  of  electrons  of 
kinetic  energy  equal  to  cut-off  energy.  Secondary  radiation  is 
also  taken  into  account  for  AD  computation.  This  secondary 
radiation  is  mainly  produced  during  the  slowing  down  of 
secondary  or  ternary  particles  ( i.e . electrons,  positrons  and 
delta  rays)  and  its  intensity  might  be  great  in  certain 
conditions  [25-26]. 

The  empiric-analytical  protocol  TG43  [1,  8-9]  is  a 
procedure  through  which  it  is  possible  to  determine  AD  rate 
in  water  around  the  seed.  By  this  formalism  and  taking  into 
consideration  the  revolution  symmetry  of  the  seed  absorbed 
dose  rate  £)  is  defined  in  spherical  coordinates  by  the 
expression: 


D(r,d ) = SkA 


GL(r,0) 

GL(r0,d0) 


gL  (r)FL  (r,  0) 


(1) 


where  SK  is  air  Kerma  strength  at  1 m distance  of  the  source 
centre.  The  dose  rate  constant  A represents  the  absorbed  dose 
rate  at  reference  distance  r0  (with  r0  = 1 cm)  on  the  mid 
transverse  plane  (0  = 0O  with  0,,  = jf2)  related  to  SK.  GL(r,  0) 
is  the  geometry  function  in  the  line  source  approximation: 


jg(r,6>)  if  0*0 
rL  sin(6>) 


GL{r,0)  = \ 1 


if  0 = 0 


(2) 


where  f is  the  angle,  in  radians,  subtended  by  the  tips  of  the 
hypothetical  line  source  with  respect  to  the  calculation  point, 
P(r,  0).  gL(r)  is  the  radial  dose  function  given  by: 

x,  _GL(ro,0o  ) D(r,  0O ) r (r„ ) D0(r) 

8lU  GL(r,0o)  D(ro,0o)  r0  D0(r0) 

where  J30(r)  = [3<r,  0„)  and  Df  r)  = D(r,0o) . Eq.  (3)  takes  into 
account  the  radial  dependence  of  absorption  and  dispersion 
of  photons  in  water  along  the  transverse  distance  from  the 
source.  The  quantity  F(r,  0)  is  the  anisotropy  function  defined 
by: 


F(r,e)=F(rA)  D(r.e)  =MLmSent0) 
GL(r,9)  D(r,6n)  P(r,0)  D0(r) 


(4) 
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Thus,  based  on  eq.  (1)  and  tabulated  data  of  A,  gL  and  F 
for  a given  seed,  it  is  possible  to  know  £)  in  the  whole  space 
(defined  by  the  phantom).  Moreover,  it  should  be  noted  from 
eq.  (3)  that  the  dependence  of  £)  with  the  r-variable  directly 
influences  the  function  gL  which  is  proportional  to  D0(r) 
multiplied  by  the  geometry  factor  r//}0(r).  Conversely,  a poor 
dependence  of  F(r,  0)  onto  the  radial  behaviour  of  £)  is 
expected.  This  statement  is  due  to  the  observed 
proportionality  between  F(r,  0)  and  dose  rate  profile 

D(r,0)/Do(r)  in  eq.  (4). 

3.  Results  and  Discussion 

In  each  simulation,  the  total  number  of  histories  simulated 
was  10n  and  the  uncertainties  achieved  were  Type  A (k=l) 
less  than  0.1  % for  distances  r<5cm  and  less  than  0.2  % 
for  distances  r < 10  cm  as  recommended  in  AAPM  and 
ESTRO  new  report  [1].  The  results  of  the  simulations  for 
both  seeds  are  presented  as  Kerma  rate  KJr)/KJr0)  and 
photons-AD  rate  profiles  D0(r)/ D0(r0)  (Fig.  3)  and  over  the 
3 material  media  (iridium,  steel  and  water).  From  the  high 
spatial  resolution  set  out  by  the  simulation,  we  are  able  to 
precisely  compute  K and  £)  on  both  sides  of  the  steel-water 
interface  (as  well  as  the  iridium-steel  interface).  AD  must 
change  continuously  when  crossing  any  interface  and  a 
detailed  analysis  of  the  steel-water  interface  is  particularly 
important  because  it  allows  us  to  accurately  know  the 
behaviour  of  AD  on  seeds  surface.  AD  continuity  in  the 
interface  is  consequence  of  the  secondary  electrons  trajectory 
not  being  cut  at  the  frontier  (secondary  electrons  deposit 
energy  on  interface  both  sides).  Differing  from  AD,  Kerma 
presents  discontinuities  on  interfaces  since  photon 
attenuation  in  each  material  medium  is  different.  Besides, 
computing  Kerma  considers  secondary  electrons  deposit 
kinetic  energy  at  the  same  place  they  are  set  in  motion. 

Charts  embedded  in  figs.  3a  ( Flexisource ) and  3b 
(; microSelectron-v2 ) show  this  different  behaviour  between 
Kerma  and  photons-AD  profiles.  We  observed  that  next  to 
the  steel  interface  in  water  phantom  ( r ^ 0.2  cm  ),  photons- 
AD  profile  oscillates  around  Kerma  profile  for  both  seeds. 
Photons-AD  profile  initially  shows  depression  with  respect  to 
Kerma  profile  very-close  to  seed  surface,  then  both  profiles 
become  equal  (in  >'  = 0.05  mm  for  both  seeds)  followed  by  an 
inverse  behaviour  in  which  photons-AD  profile  exceeds 
Kerma  profile  forming  a peak.  This  peak  in  photons-AD 
profile  over  Kerma  profile  on  seeds  vicinity  might  be 
explained  from  a reinforcement  of  the  secondary  electrons 
density  in  the  aqueous  side  of  the  steel-water  interface.  Water 
has  lower  effective  atomic  number  and  lower  density  than 
those  of  steel  and  this  causes  photons-AD  profile  to  be 
greater  than  Kerma  profile  in  the  region.  In  this  area  there  are 
backscattered  electrons  from  the  steel  (capsule)  to  the 
aqueous  medium  (phantom).  No  electronic  equilibrium  is 
established  in  the  region  nearest  the  seeds  ( r < 2 mm ). 

The  differences  that  we  found  between  Kerma  rate  and 


photons-AD  rate  profiles  for  both  seeds  are  up  to  6.4  % 
(around/'  = 0.07  mm  , see  embedded  charts  on  Figs.  3a  y 3b). 
Further  from  seeds  surface  the  differences  between  Kerma 
and  photons-AD  profiles  remained  limited  within  1 %>,  with 
Kerma  profile  below  photons-AD  profile  as  a consequence  of 
the  attenuation  of  the  primary  photons  in  water  [27]. 

The  differences  in  observed  behaviour  between  f)  and  K 
profiles  can  be  explained  by  the  lack  of  electronic 
equilibrium  near  seeds  surface,  corresponding  to  the  build-up 
region.  This  process  dominates  up  to  a distance  that  is 
equivalent  to  the  range  of  the  secondary  electrons  of  higher 
energy  ( i.e . approximately  2 mm,  for  photo-electrons  of  0.6 
MeV).  According  to  the  energy  spectrum  that  these  seeds 
emit  (Fig.  2),  the  most  energetic  photons  with  appreciable 
intensity  (see  embedded  chart  in  Fig.  2)  correspond  to  energy 
of  0.6  MeV.  These  photons  are  capable  of  producing 
secondary  electrons  with  kinetic  energy  close  to  0.6  MeV. 
Once  this  “equilibrium  thickness”  is  surpassed,  it  is  observed 
that  both  quantities  (photons-AD  and  Kerma  profiles) 
diminish  parallel  to  each  other  (see  Fig.  3)  with  photons-AD 
above  Kerma.  For  both  seeds,  the  difference  between  these 
two  quantities  (photons-AD  and  Kerma  profiles)  is  less  than 
1%  and  this  behaviour  is  observable  for  both  seeds  (as  an 
example,  see  Fig.  3c  for  Flexisource  seed).  In  this  region, 
farther  from  seeds  surface  the  small  differences  observed 
between  Kerma  and  photons-AD  profiles  indicate  that  the 
primary-photons  attenuation  is  a process  that  could  be 
compensated  partially  by  another  one,  originating  in  the 
secondary-electrons  bremsstrahlung-  losses.  Primary -photons 
attenuation  helps  photons-AD  to  be  greater  than  Kerma, 
photons-AD  at  one  point  would  equal  Kerma  at  another  point 
at  which  primary-photons  density  is  greater.  On  the  other 
hand,  radiation  losses  diminish  photons-AD  but  not  Kerma. 
Flence,  primary-  photons  attenuation  process  would  be  barely 
more  important  than  secondary-particles  radiation,  making 
photons-AD  profile  to  prevail  over  Kerma  profile,  although 
just  in  a small  amount. 

As  it  has  been  mentioned  before  (see  introduction  section) 
several  authors  have  studied  build-up  region  in  the  past  for 
FIDR-brachytherapy  seeds.  These  computations  were  based 
on  simplified  seeds  modelling,  either  in  seed  design  and/or  in 
the  energy  spectrum  emitted  by  it.  Ballester  et  al  [4] 
calculated  AD  in  water  from  spherical  Ir-192  seeds  in  a steel 
capsule.  These  spheres  have  an  active  core  diameter  and 
capsule  thickness  comparable  to  transversal  seed  dimensions. 
They  found  an  oscillating  behaviour  between  photons-AD 
and  Kerma  similar  to  the  one  found  in  this  work,  where  no 
approximations  have  been  made.  Luxton  and  Joszef  [6] 
evaluated  AD  in  water  by  mono-energetic  photon-emitting 
point  sources  ranging  from  10  keV  to  2 MeV.  They  found 
build-up  regions  of  1 mm  for  500  keV  sources  and  2 mm  for 
800  keV  sources.  This  result  is  also  consistent  with  the  one 
achieved  for  the  seeds  simulated  in  this  work  without 
simplifications,  either  in  seeds  geometry  or  in  seed  photon- 
emitting  spectrum. 

Table  II  displays  the  results  achieved  in  this  work  for  both 
seeds  Kerma  values  relative  to  photons-AD.  Results  achieved 
in  this  work  for  microSelectron-v2  seed  are  compared  to 
those  achieved  by  other  authors  in  the  past  [4,  11,  28].  On  the 
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right  of  each  column  that  corresponds  to  quotient 
K(r)/i)(r)  values  achieved  by  other  authors,  another 
column  with  percentage  differences  8 % between  our  results 
and  those  from  these  authors  is  added.  From  table  II,  general 
agreement  among  different  authors  results  may  be  observed. 
The  differences  found  between  this  work  results  and  those  of 
other  authors  are  mostly  less  than  or  equal  2 %.  For 
Flexisource  seed,  it  has  not  been  possible  to  compare  our 
results  of  kerma  rate  to  dose  rate  values  because  of  the  lack 
of  other  authors  results  of  photons-AD  and  Kerma  in  the 
near-seed  region. 

From  eq.  (3)  and  K and  [)  profiles  we  computed  radial 
dose  functions  from  Kerma  ( g,  K0~)  ) and  from  photons-  AD 
( Sl.d(t ) ) for  both  seeds  (see  Fig.  4).  For  both  seeds,  it  is 
observed  that  gL,D(r)  clearly  deviates  from  gL.K(r)  up  to 


around  2 mm  distance.  This  distance  corresponds  to  the 
equilibrium  thickness  (see  embedded  charts  in  Figs.  4a  and  4b). 
The  differences  between  gL.K(r)  and  gLiD(r)  in  this  region,  up 
to  6.4  % are  a consequence  of  the  same  relative  differences  on 
photons-AD  rate  and  Kerma  rate  profiles.  The  function 
obtained  from  our  data  continues  very  close  to  the  one 
obtained  from  photons-AD  rate  profile  although  gLD(r)  values 
stay  just  above,  this  difference  being  less  than  or  equal  1 %>. 
Results  achieved  in  this  work  for  g(r)  are  consistent  with  the 
ones  calculated  in  the  past  by  Granero  et  al  [2],  Daskalov  et  al 
[3]  corrected  to  full  scattering  conditions  [1,  29]  and  Taylor 
and  Rogers  [11].  Table  III  in  the  appendix  contains  values 

achieved  in  this  work  for  dose  rate  profile  D0(r)l  D0(r0)  and 
gL,D(r)  photons-AD  radial  dose  function  and  are  presented 
with  spatial  resolution  of  0.05  mm  up  to  5 mm. 


(c)  io'3  1CT2 


103 


10 3 r (cm)  10 2 
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Figure  3.  Dose  rate  profiles  for  6 = 6q  calculated  from  photons-AD  rate  (black  scatter)  and  Kerma  rate  (gray  scatter):  (a)  Flexisource  results  from  this  work 
and  (b)  microSelectron-v2  results  from  this  work.  Embedded  window  chart  shows  amplification  of  photons-AD  and  Kerma  rate  profiles  near  steel-water 
interface  in  figures  (a)  and  (b).  In  (c),  amplification  shows  photons-AD  rate  and  Kerma  rate  profiles  in  three  different  segments:  inside  the  seed  (iridium  core 
and  steel  capsule),  near  seed  capsule  and  water  phantom  for  Flexisource  seed. 
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Table  II.  Kerma  values  related  to  photons-AD  achieved  in  this  work  are  compared  to  those  obtained  by  other  authors  in  the  past:  Wang  and  Li  [27],  Taylor 
and  Rogers  [10],  Ballester  et  al  [4]  for  microSelectron-v2  seed.  On  the  right  of  each  column  of  results  from  other  authors,  another  8 % percentage  difference 
column  (between  this  work  and  other  author  values)  is  added.  Only  this  work  results  are  presented  for  Flexisource  seed  because  of  the  lack  of  Kerma  and 
photons-AD  rate  near-source  results  from  other  authors. 


r (mm) 

K(r)/D(r) 

Flexisource 

microSelectron  v2 

This  work 

Wang  and  Li  2002 

8 (%) 

Taylor  and  Rogers 
2008 

8 (%) 

Ballester  et  al  _ 

2009  5(/o) 

This  work 

0.44 

1.059 

- 

- 

- 

- 

- 

** 

- 

0.45 

1.045 

- 

- 

- 

- 

- 

- 

- 

0.46 

1.033 

- 

- 

- 

- 

- 

- 

1.018 

0.47 

1.022 

- 

- 

- 

- 

- 

- 

1.008 

0.48 

1.013 

- 

- 

- 

- 

- 

- 

0.999 

0.49 

1.004 

- 

- 

- 

- 

0.992 

0.50 

0.997 

1.035 

5.2 

- 

- 

0.997 

1.3 

0.984 

0.55 

0.970 

- 

- 

- 

- 

0.958 

0.60 

0.957 

0.961 

2.0 

0.958 

1.7 

0.973 

3.3 

0.942 

0.65 

0.951 

- 

- 

- 

- 

- 

- 

0.934 

0.70 

0.949 

0.942 

1.2 

0.941 

1.1 

0.941 

1.1 

0.931 

0.75 

0.949 

» 

- 

- 

- 

- 

- 

0.932 

0.80 

0.952 

0.940 

0.6 

0.943 

1.0 

0.935 

0.1 

0.934 

0.85 

0.955 

- 

- 

- 

- 

- 

- 

0.938 

0.90 

0.959 

0.946 

0.4 

0.952 

1.1 

0.940 

0.2 

0.942 

0.95 

0.964 

--- 

- 

- 

- 

- 

- 

0.946 

1.00 

0.968 

0.959 

0.9 

0.960 

1.1 

0.951 

0.1 

0.950 

1.10 

0.976 

- 

- 

0.961 

0.3 

0.958 

1.20 

0.982 

- 

- 

0.978 

1.6 

0.972 

0.9 

0.963 

1.30 

0.986 

- 

- 

0.978 

1.0 

0.968 

1.40 

0.989 

- 

- 

0.986 

1.5 

0.982 

1.1 

0.971 

1.50 

0.990 

- 

- 

0.986 

1.1 

0.975 

1.60 

0.991 

- 

- 

0.990 

1.3 

0.988 

0.9 

0.977 

1.70 

0.992 

- 

- 

0.992 

1.2 

0.980 

1.80 

0.991 

- 

- 

0.995 

1.3 

0.995 

1.3 

0.982 

1.90 

0.992 

- 

0.995 

1.1 

0.984 

2.00 

0.994 

0.991 

0.5 

0.999 

1.3 

0.994 

0.8 

0.986 

4.  Conclusions 

For  this  work  we  have  selected  two  of  the  most  commonly 
used  seeds,  which  remain  valid  today  in  many  medical 
centres  for  HDR-brachytherapy  treatments.  The  selected 
seeds  (Flexisource  and  microSelectron-v2 ) come  from 
different  manufacturers,  giving  greater  generality  to  the 


conclusions  achieved  in  this  work. 

On  the  other  hand  it  should  also  be  noted  that  seeds 
modelling  and  cross  sections  used  in  our  simulations  were 
the  same  for  both:  Kerma  rate  and  photons-AD  rate 
calculations.  Thus  the  observed  differences  in  behaviour 
between  Kerma  and  photons-AD  profdes  are  attributed 
exclusively  to  partial  loss  of  electronic  equilibrium  in 
different  regions  of  the  space  surrounding  the  seeds. 
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Figure  4.  Radial  dose  function  calculated  from  photons-AD  rate  gL.o(r)  (black  line)  and  Kerma  rate  gL,tc(r)  (gray  line)  are  presented.  In  (a),  this  work 
Flexisource  data  is  compared  to  the  ones  obtained  by  Granero  et  al  [2]  in  green  diamonds.  In  (b),  this  work  microSelecetron-v2  data  is  compared  to  those 
obtained  by  Daskalov  et  al  [3]  in  blue  squares1"1,  Granero  et  al  [2]  in  green  diamonds  Taylor  and  Rogers  [11]  in  red  diamonds  ,**),  and  Granero  et  al  [12] 
in  black  line  and  white  filled  circles  r \ (*)  Achieved  from  Kerma  rate  profile,  and  corrected  to  full  scattering  conditions  [1,  29],  (f)  From  Kerma  rate  profile. 
(**)  From  AD  rate  profile,  (ff)  Mixed  form  AD  rate  and  Kerma  rate  profiles,  and  for  microSelectron-v2r  seed.  Granero  et  al  [12]  AD  profile  dose  rate 
function  includes  and  electrons  emitted  from  source. 


Phantom  discretization  into  small  voxels  (see  Section  2) 
was  due  to  the  need  for  greater  spatial  resolution  than  that 
achieved  in  the  past  according  to  the  literature  of  the  subject, 
mainly  in  the  area  of  high  spatial-gradient  close  to  the  seed. 
We  found,  as  expected,  a continuous  behaviour  of  photons- 
AD  across  the  steel-water  interface  opposite  to  Kerma  which 
showed  a discontinuity  in  the  same  interface.  By  using  high 
spatial  resolution  and  the  continuity  of  AD  through  the 


interface,  it  was  possible  to  calculate  the  precise  value  of 
photons-AD  on  seeds  surface.  Analysing  these  results  on  the 
three  material  media  ( iridium , steel  and  water)  and 
particularly  in  the  steel-water  interface,  we  found  that  we 
were  able  to  calculate  AD  with  precision  in  the  build-up 
region.  Whenever  comparisons  are  possible,  our  results  are  in 
agreement  with  those  reported  by  other  authors. 

Quantities  gL,  F and  A completely  characterize  delivered 
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dose  in  water  from  Iridium-192  seeds  that  are  commonly 
used  in  HDR-brachytherapy,  since  they  allow  obtaining  AD 
in  water  in  the  space  surrounding  the  seeds.  From  those 
quantities  (by  its  definition  in  TG43  formalism),  function 
gL(f)  can  be  seen  to  be  directly  affected  by  radial  behaviour 
of  absorbed  dose  profile. 
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Appendix 


Table  III.  Both  seeds  photons-AD  rate  profile  D0(r)/D0(r0)  (in  a)  and  radial  dose  function  gL,D(r)  calculated  from  photons-AD  rate  profile  (in  b)  achieved  in 
this  work  are  presented  in  r-steps  of  0.05  mm  up  to  a distance  of  5 mm.  Each  column  contains  data  for  both  simulated  seeds:  Flexisource  (on  the  left)  and 
microSelectron-v2  (on  the  right). 


a) 


i>o(r)/DM 

r (mm) 

0 

1 

2 

3 

4 

.00 

60.748 

61.457 

20.423 

20.595 

10.041 

10.097 

5.9015 

5.9202 

.05 

56.522 

57.122 

19.587 

19.752 

9.7445 

9.7975 

5.7651 

5.7830 

.10 

52.737 

53.241 

18.800 

18.959 

9.4610 

9.5106 

5.6333 

5.6504 

.15 

49.332 

49.755 

18.059 

18.211 

9.1893 

9.2359 

5.5059 

5.5223 

.20 

46.257 

46.615 

17.361 

17.506 

8.9290 

8.9727 

5.3828 

5.3984 

.25 

43.470 

43.777 

16.701 

16.840 

8.6792 

8.7203 

5.2637 

5.2786 

.30 

40.934 

41.202 

16.077 

16.210 

8.4396 

8.4783 

5.1484 

5.1627 

.35 

38.616 

38.858 

15.486 

15.614 

8.2096 

8.2460 

5.0368 

5.0505 

.40 

36.491 

36.715 

14.928 

15.049 

7.9887 

8.0230 

4.9288 

4.9419 

.45 

154.82 

162.78 

34.536 

34.751 

14.399 

14.514 

7.7764 

7.8087 

4.8241 

4.8368 

.50 

142.75 

147.84 

32.734 

32.943 

13.897 

14.006 

7.5723 

7.6027 

4.7228 

4.7349 

.55 

130.79 

134.24 

31.068 

31.273 

13.421 

13.523 

7.3759 

7.4047 

4.6245 

4.6361 

.60 

119.48 

121.99 

29.523 

29.728 

12.968 

13.065 

7.1869 

7.2141 

4.5293 

4.5404 

.65 

109.04 

111.02 

28.090 

28.292 

12.538 

12.629 

7.0049 

7.0307 

4.4369 

4.4476 

.70 

99.575 

101.23 

26.757 

26.952 

12.129 

12.214 

6.8295 

6.8541 

4.3473 

4.3575 

.75 

91.086 

92.512 

25.515 

25.709 

11.739 

11.818 

6.6605 

6.6840 

4.2603 

4.2702 

.80 

83.512 

84.764 

24.356 

24.548 

11.367 

11.441 

6.4976 

6.5200 

4.1759 

4.1854 

.85 

76.778 

77.878 

23.274 

23.462 

11.012 

11.082 

6.3405 

6.3619 

4.0940 

4.1031 

.90 

70.797 

71.758 

22.261 

22.445 

10.674 

10.739 

6.1593 

6.1796 

4.0144 

4.0231 

.95 

65.482 

66.311 

21.313 

21.491 

10.350 

10.411 

6.0427 

6.0622 

3.9371 

3.9455 

b) 


gL.p(r) 


r (mm) 

0 

1 

2 

3 

4 

.00 

1.0124 

1.0171 

0.9927 

0.9927 

0.9938 

0.9927 

0.9948 

0.9942 

.05 

1.0085 

1.0136 

0.9928 

0.9924 

0.9938 

0.9928 

0.9948 

0.9943 

.10 

1.0049 

1.0106 

0.9928 

0.9922 

0.9939 

0.9929 

0.9949 

0.9943 

.15 

1.0019 

1.0081 

0.9929 

0.9921 

0.9939 

0.9930 

0.9949 

0.9944 

.20 

0.9994 

1.0059 

0.9929 

0.9919 

0.9940 

0.9931 

0.9950 

0.9944 

.25 

0.9973 

1.0041 

0.9930 

0.9918 

0.9940 

0.9932 

0.9950 

0.9945 

.30 

0.9958 

1.0025 

0.9930 

0.9917 

0.9941 

0.9933 

0.9951 

0.9945 

.35 

0.9946 

1.0012 

0.9931 

0.9916 

0.9941 

0.9934 

0.9951 

0.9946 

.40 

0.9938 

1.0001 

0.9932 

0.9915 

0.9942 

0.9935 

0.9952 

0.9946 

.45 

0.9621 

0.9358 

0.9933 

0.9991 

0.9932 

0.9915 

0.9942 

0.9936 

0.9952 

0.9947 

.50 

0.9908 

0.9779 

0.9929 

0.9982 

0.9933 

0.9915 

0.9943 

0.9937 

0.9953 

0.9947 

.55 

1.0102 

1.0054 

0.9927 

0.9973 

0.9933 

0.9916 

0.9943 

0.9937 

0.9953 

0.9948 

.60 

1.0222 

1.0221 

0.9926 

0.9966 

0.9934 

0.9916 

0.9944 

0.9938 

0.9954 

0.9948 

.65 

1.0288 

1.0310 

0.9926 

0.9959 

0.9934 

0.9917 

0.9944 

0.9939 

0.9954 

0.9949 

.70 

1.0314 

1.0345 

0.9924 

0.9952 

0.9935 

0.9919 

0.9945 

0.9939 

0.9954 

0.9949 

.75 

1.0310 

1.0346 

0.9924 

0.9946 

0.9935 

0.9920 

0.9945 

0.9940 

0.9955 

0.9950 

.80 

1.0287 

1.0323 

0.9925 

0.9941 

0.9936 

0.9921 

0.9946 

0.9940 

0.9955 

0.9950 

.85 

1.0252 

1.0289 

0.9925 

0.9937 

0.9936 

0.9922 

0.9946 

0.9941 

0.9956 

0.9951 

.90 

1.0211 

1.0249 

0.9926 

0.9933 

0.9937 

0.9924 

0.9947 

0.9941 

0.9956 

0.9951 

.95 

1.0167 

1.0209 

0.9927 

0.9929 

0.9937 

0.9925 

0.9947 

0.9942 

0.9957 

0.9952 
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